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MODE I, MODE II, AND MIXED-MODE FRACTURE OF 
PLASMA-SPRAYED THERMAL BARRIER COATINGS 
AT AMBIENT AND ELEVATED TEMPERATURES 


Sung R. Choi 
Ohio Aerospace Institute 
Brook Park, Ohio 44142 

Dongming Zhu 

U.S. Army Research Laboratory 
Glenn Research Center 
Cleveland, Ohio 44135 

Robert A. Miller 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

The mixed-mode fracture behavior of plasma-sprayed Zr 02-8 wt% Y 2 O 3 thermal 
barrier coatings was detennined in air at 25 and 1316 °C in asymmetric four-point flexure 
with single edge v-notched beam (SEVNB) test specimens. The mode I fracture 
toughness was found to be A! Ic = 1.15 ± 0.07 and 0.98 ± 0.13 MPa 4m , respectively, at 25 
and 1316 °C. The respective mode II fracture toughness values were Al llc = 0.73 ± 0.10 
and 0.65 ± 0.04 MPaVm . Hence, there was an insignificant difference in either A! Ic or 
A! IIc between 25 and 1316 °C for the coating material, whereas there was a noticeable 
distinction between K lc and K Uc , resulting in K Uc /Ki c = 0.65 at both temperatures. The 
empirical mixed-mode fracture criterion best described the coatings’ mixed-mode 
fracture behavior among the four mixed-mode fracture theories considered. The angle of 
crack propagation was in reasonable agreement with the minimum strain energy density 
criterion. The effect of the directionality of the coating material in on K ic was observed to 
be insignificant, while its sintering effect at 1316 °C on K ic was significant. 


1. INTRODUCTION 

Thermal barrier coatings (TBCs) have attracted increasing attention for advanced gas 
turbine and diesel engine applications due to their ability to provide thermal insulation to 
engine components . 1-3 The merits of using the ceramic thermal barrier coatings are well 
recognized and include the potential increase in engine operating temperature with 
reduced cooling requirements, resulting in significant improvement in thermal efficiency, 
performance, and reliability. Plasma-sprayed zirconia-based ceramics are one of the most 
important coating materials in light of their low thermal conductivity, relatively high 
thermal expansivity, and unique microstructure as a result of the plasma spraying process. 
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However, the limited durability of thermal barrier coatings under severe thermal and 
mechanical loading conditions encountered in heat engines remains one of the major 
problems. As a result, the development of thermal barrier coatings requires a better 
understanding of both thermal and mechanical behavior of the coating materials to ensure 
life and reliability of the related components. 

It has been suggested that the important limiting factor encountered in thin plasma- 
sprayed thermal barrier coatings is the relatively low fracture energy of the coating in 
planes close to and parallel to the interface. Various efforts have been made to determine 
mode I ‘interfacial’ fracture toughness of the coatings in the vicinity of the interface 
using various techniques such as the indentation method, debonding technique, three- or 
four-point flexure delamination technique, and compact tension test method . 4-6 Failure of 
the thick thermal barrier coatings has been observed within the bulk of the coating 
material, independent of any delamination that typically occurs at the interface of a thin 
coating and a substrate . 7 There have been efforts to determine mode I fracture toughness 
of thick thermal barrier coatings at ambient and elevated temperatures. In fact, the 
majority of data on fracture toughness of both thin and thick coating materials have been 
determined under mode I loading. However, rarely are structural components or coatings 
subject to pure mode I loading. This is particularly true for thermal barrier coatings that 
encounter complex thermal and mechanical loading in engine operations. Recently, 
Callus and Berndt 8 used a pure-shear technique to determine the interfacial critical mode 
II strain energy release rate of some thin coatings. Their data, however, were only for 
ambient temperature. Mixed-mode data on either thin or thick coatings are rarely 
available in the literature at elevated temperatures, despite an important fact that 
mechanical behavior of coating materials should not be based solely on ambient- 
temperature properties. Dense oxide ceramics (such as zirconia) have exhibited 
significant changes in mechanical properties such as strength, fracture toughness, and 
slow crack growth when temperature is increased above 800 to 1000 °C . 9 

The objective of this work was to determine modes I and II fracture toughness and 
mixed-mode fracture behavior of free-standing thermal barrier coatings of plasma- 
sprayed Z 1 O 2-8 wt% Y 2 O 3 at an elevated temperature of 1316 °C in air. The choice of 
this temperature was based on a typical target temperature of aerospace gas turbine 
applications. These same properties were also determined at ambient temperature 
(25 °C). An asymmetric four-point flexure test technique was used at both temperatures 
in conjunction with single-edge-v-notched beam (SEVNB) test specimens, which yielded 
simplicity in both specimen and crack preparation and in test procedure. Fracture locus, 
ranging from mode I to mode II, was determined. Several mixed-mode fracture criteria 
were analyzed based on experimental data on mode I and mode II stress intensity factors 
as well as on crack propagation angles. The mixed-mode fracture behavior of the coatings 
was compared with that of typical advanced structural ceramics in terms of mixed-mode 
fracture theories. The effects of sintering and the material’s directionality on mode I 
fracture toughness at ambient temperature were also characterized. 
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2. EXPERIMENTAL PROCEDURES 


2.1. Material 

The Zr02-8 wt% Y2O3 powder with an average particle size of 60 pm was first 
plasma-sprayed on a graphite substrate measuring 150 by 100 by 6.5 mm to a thickness 
of about 6 mm, using a Sulzer-Metco ATC-1 plasma coating system with an industrial 
robot. The plasma-spray conditions can be found elsewhere. 10 A free standing, plasma- 
sprayed ceramic billet was then obtained by burning away the graphite substrate at 
680 °C in air for 24 h. The billet was machined into the final, rectangular flexure test 
specimen with nominal dimensions of 3 by 4 by 50 mm, respectively, in width, depth, 
and length. The 3-mm-wide face of flexure test specimens corresponded to the plane 
perpendicular to the plasma spraying direction. Major physical and mechanical properties 
of as-sprayed coating material including hardness, fracture toughness, and strength have 
been determined previously at ambient temperature 11 ' 12 and are presented in Table 1. 
Figure 1 shows a typical fracture surface and a polished surface showing the 
microstructure of as-processed coatings, in which large amounts of microcracks and 
pores are characterized in conjunction with a unique platelet structure. 

2.2. Mixed-Mode Fracture Testing 

2.2.1. Preparation of sharp precracks 

Sharp v notches were introduced in flexure test specimens, using the single-edge-v- 
notched-beam method. 13 This method utilizes a razor blade with diamond paste to 
introduce a sharp root radius by tapering a saw cut. Sharp v-notch radii ranging from 4 to 
6 pm have been successfully obtained for alumina, glass ceramic, silicon nitride, 
zirconia, and silicon carbide ceramics. 13 A starter straight-through notch 0.6 mm deep and 
0.026 mm wide was made on the 3-mm-wide face of the test specimens. A steel razor 
blade was put into the starter notch sprinkled with diamond paste with a particle size of 
9 pm. Typically, a load of about 10N was applied through the razor blade with a stroke 


Table 1. Typical Physical and Mechanical Properties of Plasma-Sprayed Zro2-8 Wt% 
Y2O3 Thermal Barrier Coatings at Ambient Temperature 11 ' 12 


Density, 

g/cm 3 

Elastic modulus 3 
(in compression), 
E, 

GPa 

Fracture 
toughness, b 

Strength 


K lc , 

MPa Vrn 

Type of test 

Number 

of 

specimens 

Average 
strength, b 
MPa 

5.22 

34 

1 .0(0.1) 

Tension 

10 

10 (2) 




Trans-thickness 

10 

11(1) 




tension 






compression 

10 

324(72) 




Four-point flexure 

20 

32(6) 




Biaxial flexure 

10 

40(4) 


■“The value of E is estimated from the starting point of the stress-strain curve. 1 1 
b The number in parenthesis represents the +1.0 standard deviation. 
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(a) 50 nm (b ) 10 nm 

Figure 1. — Typical microstructures of plasma-sprayed Zr0 2 -8 wt% Y 2 0 3 thermal barrier coatings, 
(a) Polished surface, (b) Fracture surface. Arrow indicates plasma-spray direction. 


rate (back-and-forth motion) of 1 Hz and a stroke distance of about 13 mm in a specially 
designed polishing machine. 11 The final notch depth and root radius were 2.0 mm and 20 
to 50 |im, respectively, resulting in the crack size-to-specimen depth ratio of 0.5 (=a/W). 
The coating material resulted in less sharpness in root radius due to its porous and 
microcracked nature, compared to typical dense ceramics. However, it has observed that 
the sharpness ranging from 20 to 50 pm was sufficient to give a consistent and accurate 
value of fracture toughness of the coating material. 11 Note that the through-the-thickness 
sharp notches thus prepared were aligned parallel with respect to the plasma spraying 
direction. 

Other methods to generate sharp cracks to estimate fracture toughness, such as the 
single-edge-precracked-beam (SEPB) 14 and the indentation techniques, were not feasible 
for the coating material: The indentation response was poor because of the material’s 
significant porosity, microcracks, and ‘softness,’ so that well-defined indentation cracks 
were not achieved. This indicates that microindentation techniques may not be 
appropriate for plasma-sprayed coating materials whose pore and microcrack sizes in 
some cases are much greater than microindentation cracks, and as a result, a continuum 
approach is no longer valid. 11 Also note that cracks produced by this SEVNB method 
were very uniform in their crack-front geometry and very consistent in size with high 
reproducibility from specimen to specimen. This gives rise to very little scatter in or a 
small coefficient of variation in the values of fracture toughness, a notable advantage of 
the SEVNB method over other methods. 

2.2.2. Asymmetric four-point flexure testing 

The sharp v-notched flexure test specimen was loaded in asymmetric four-point 
flexure as shown in figure 2. The ratio of mode I to mode II stress intensity factor (SIF), 
Ki/Kn, was varied by varying the distance of the precrack from the center plane, s, as 
shown in the figure. When the precrack was centered with respect to the loading point 
(s = 0), the precrack was subjected to pure mode II loading. As .v was increased, K\/K u 
also increased. The stress intensity factors, K\ and K n , are expressed as 
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Figure 2. — Schematic of asymmetric four-point flexure 
geometry with accompanying shear force and bending 
moment diagrams. 


K^aina) 112 ^^ ( 1 ) 

K u =x{na) m F ll [^ ( 2 ) 

where <7 is the applied (remote) normal stress, T is the applied (remote) shear stress, a is 
the crack size, and W is the specimen depth. F\ and F u are crack geometry factors in 
modes I and II, respectively. The normal and shear stresses are given from the elementary 
beam theory by 


A-B ( 6 sP \ 
A+B\bW 2 ) 


( 3 ) 


A-B 

x = 

A + B 



(4) 


where A and B are distances from a load point as defined in figure 2, b is the specimen 
width, and P is the applied force. The values of A = 10 mm and B = 5 mm were typically 
used, but in some cases ( Kn c at RT) values of A = 12 and B = 6 mm were used for 
comparison and verification. The test specimen’s cross section was W = 4 mm and b = 
3 mm. 
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— o — F|j Suresh [15] 
— * — Fn; Suresh [15] 
F|| Wang [16] 

— a — F||; Wang [16] 

F|| He[1 7] 

F,i; He [17] 

F|| Srawley [19] 


Asymmetric four-point flexure 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


a/W 


Figure 3. — Various solutions of crack geometry factors F| and F|| applied in 
the asymmetric four-point flexure configuration. 


Several different expressions of F: and F lr were suggested by Suresh et al., Wang 
etal., 16 and Fie and Hutchinson 17 for the case of through-the-thickness crack, and their 
respective results of F : and F n are illustrated in figure 3. As seen from the figure, no 
significant difference in F: and F n between the solutions was found, particularly when 
a/W is between 0.35 and 0.50. The solution by He and Hutchinson 17 that provided a 
convenient polynomial expression was used in here, and is given as follows: 


f {w) = 1122 - l, 21 (t) + 3 - 74 <t) 2 + 3 ' 873 (w ) 3 


_( a 'i 4 ( 

a V 

) — + 22.55 

— 

Vw) l 

.WJ 


w 


( 5 ) 


fJ — ) = 7.264 - 9.37 ( — ) + 2.74f — 


+ 1.87f— 1 1 for 0 < — < 1 

VWJ KWJ W 


(6) 


Note that Eq. (5) was quoted from Murakami’s ls and was almost identical to Srawly and 
Gross’ solution 19 up to a/W= 0.7, as seen in the figure. Also note that the difference in F: 
and Fn between solutions in the range of a/W= 0.4 to 0.5 was negligible, resulting in 
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F \ = /’n 

Hence, in the case of a/W = 0.4 to 05, the K u /Ki ratio (or called ‘mixity’) can be 
simplified from Eqs. (1)— (4) and (7) as follows: 

K R= WFu_^W 
Ki 6sF l 6s 


(3 = tan 1 


f 


V 


Ku 


K, 


I ) 



(9) 


where (3 is the mixity angle. The angle (3 was also called the “equivalent crack angle” by 
Maccagno and Knott 20 for situations which do not actually use the inclined crack 
configuration. The ratio of crack size to specimen width used in this work was chosen to 
be a/W = 0.5 so that Eq. (8) or (9) could be conveniently used. The values of 5 = 0 to 
3.6 mm were selected to give uniformly distributed (3 between pure mode I ((3 = 0°) and 
pure mode II ((3 = 90°) based on Eq. (9). Pure mode I fracture toughness K Ic was 
determined using a symmetric four-point flexure fixture with 10-mm inner and 20-mm 
outer distances. 

All testing was performed in displacement control using a SiC flexure fixture with an 
actuator speed of 0.5 mm/min of an Instron electromechanical test frame (Model 8562, 
Instron, Canton, MA). Testing was conducted at 25 and 1316 °C in air. Each test 
specimen in elevated-temperature testing was held for thermal equilibration for about 
20 min prior to testing. At room temperature, four and nine specimens were tested for K\ c 
and K llc , respectively, whereas at 1316 °C, four specimens were examined for each test. 
There were nine mixity values associated with the specimens tested at 25 °C, and four 
mixity values associated with those tested at 1316 °C. Typically, two or three specimens 
of each mixity value were tested. After testing, the crack size of each tested specimen 
was determined optically from its fracture surface based on the three-point measurements 
in accordance with test standard ASTM Cl 421. 14 The crack propagation angles were also 
determined. 


2.3. Effects of Sintering and Material Directionality 

In order to determine the effect of sintering on mode I fracture toughness of the 
coating material, additional fracture toughness testing was carried out at ambient 
temperature using the test specimens annealed at 1316 °C in air for 5, 20, 100, and 500 h. 
The sharp v-notches were introduced in the annealed specimens using the same technique 
described in section 2.2.1. The test fixture, test frame, and test procedure used were the 
same as those used in the aforementioned ‘regular’ fracture toughness testing. Fracture 
surfaces were also examined via SEM with respect to annealing time to see any change in 
microstructure and morphology. 

Mode I fracture toughness (K Ic ) of as-sprayed coating material was also evaluated at 
ambient temperature in the plane perpendicular to the plasma-sprayed direction (that is, 
within the ‘layered’ direction) using the double cantilever beam (DCB) method. A sharp 
v notch was introduced in the center plane of the 4-mm-wide face of each flexure test 
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◄ — 

► K\ 

=//= 

==/ 


(b) 


Figure 4. — (a) A schematic geometry of double-cantilever-beam (DCB) test specimens used to 
determine K\ c in the plane perpendicular to the plasma-sprayed (PS) direction, (b) ‘Regular’ 
fracture toughness testing by the SEVNB method with crack planes parallel to the plasma- 
spraying direction is shown for comparison. 


specimen with a crack length of about 6 mm (see fig. 4). Fracture toughness was 
calculated based on the formula by Murakami. 18 A total of three DCB specimens were 
tested. This additional fracture testing was to determine directionality of the material in 
response to fracture toughness relative to the plasma-spraying direction. 


3. RESULTS AND DISCUSSION 
3.1. Test Results 

The results of the mode I, mode II, and mixed-mode fracture testing for the coating 
material at both 25 and 1316 °C are presented as K n versus K\ in figure 5. Each point 
represents a single datum, and the values K ic and K nc represent the average of four and 
eight measurements, respectively, at 25 °C and the average of four measurements for 
each at 1316 °C. The effect of test distance on A IIc at 25 °C was found to be insignificant: 
Kn c = 0.74±0.1 and 0.69±0.09 MPaVm, respectively, for AtB values of 10/5 and 
12/6 mm (see fig. 2). Hence, the 10/5-mm distances were exclusively used throughout the 
mixed-mode test program. At 25 °C, K\ c = 1.15±0.07 and A IIc = 0.73+0.10 MPa Vm , and 
at 1316 °C, K\ c = 0.98±0.13 and A IIc = 0.65±0.04 MPaVm . Note that the values of A IIc 
were 37 and 34 percent lower than those of K\ c at 25 and 1316 °C, respectively. It has 
been reported that for some dense ceramics A IIc was greater than A/ Ic , presumably 
attributed to the frictional interaction between the two crack planes. 21 However, the 
previous studies on advanced (dense) ceramics including silicon nitrides, alumina and 
zirconia 2 - 23 showed a different result that K nc was almost identical to K lc for a given 
material, indicative of an insignificant frictional effect on A IIc by either coarse grained or 
fine grained ceramics. The coating material, however, did not exhibit a similar value in 
both A IIc and A Ic but rather yielded a lower value in A IIc than in K lc at both 25 and 
1316 °C. It is believed that because of the porous, microcracked, ‘soft’ nature the coating 
material, there would not have been any significant frictional contribution to K Uc - 
Furthermore, the coating material would be easier to shear (mode II) than to cleave 
(mode I), probably because of its unique platelet structure formed by plasma-spraying 
process, which results in low K lic compared with A Ic . The difference between A IIc and A Ic 
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Figure 5. — A summary of K\\ versus Kj for mode I, mode II, and mixed-mode 
fracture of plasma-sprayed ZrC> 2-8 wt% Y 2 O 3 thermal barrier coatings at 
both 25 and 1316 °C in air. 


for various ceramics and brittle materials will be described in more detail in a next 
section. 

As shown in figure 5, because K nc < A) c , the overall versus K\ relation is 
represented by an ellipse at either 25 or 1316 °C with K lc being a major axis. For a given 
mixity angle ranging from (3 = 0 to 90°, K E q [= (/f [ 2 +A II 2 ) 12 ] was lower at 1316 °C than at 
25 °C, probably because of a somewhat softening effect at elevated temperature. The 
values of K u and K Uc were decreased by 15 and 11 percent, respectively, when 
temperature was increased from 25 to 1316 °C. Flowever, this decrease was marginal, as 
can be seen in a summary of as a function of temperature in figure 6, in which the Aj c 
data previously determined 1 1 at 25 and 800 °C by the SEVNB method were included with 
the current data. K\ c remains statistically invariable with an average of K\ c = 
1.0±0.1 MPaVm , regardless of test temperature up to 1316 °C. The same would be true 
for K Ec (=0.69 MPaVm) as depicted in the figure, notwithstanding the lack of 
intermediate-temperature data. The insensitivity of K\ c to temperature was also observed 
in a similar but dense (hot-pressed) 10 mol% yttria-stabilized zirconia (10-YSZ) in which 
fracture toughness of the material — evaluated in air by the same SEVNB method — was 
almost unchanged {K\ c = 1. 6-1. 8 MPaVm ) from 25 to 1000 °C. 9 

Figure 7 shows typical examples of crack propagation with different mixities ranging 
from pure mode II ( K\/Kn = 0) to pure mode I (K\IK U = 00 ). The angle of crack 
propagation 0 C was defined as an angle of the crack plane propagated with respect to the 
original precrack plane. The direction of crack plane propagated was in most cases 
straight with some exceptions such as tortuous or not-straight paths. The angle of crack 
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Figure 6. — Fracture toughnesses of K\ c and K\\ c as a function of temperature 
for plasma-sprayed Zr02-8 wt% Y2O3 thermal barrier coatings. The previous 
K\ c data 11 were included. Error bars represent ±1 .0 standard deviation. 
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(a) (b) 

1 
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(c) (d) 

Figure 7. — Typical examples of specimens showing angles of crack prop- 
agation with respect to K\/K\\ of plasma-sprayed Zr02-8 wt% Y2O3 
thermal barrier coatings: (a) pure mode II; K\!K\\ = 0, (b) K\/K\\ = 2.7, 

(c) K|/K|| = 5.5, and (d) pure mode I; K\!K\\ = =0. 
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Figure 8. — Crack propagation angles as function of K\!K\\ of plasma- 
sprayed Zr0 2 -8 wt% Y 2 0 3 thermal barrier coatings determined in 
asymmetric four-point flexure at both 25 and 1 31 6 °C in air. 


propagation was greatest in pure mode II, and decreased with increasing K { !K U reaching 
to 9 C = 0 in pure mode I. A summary of experimental data on the angle of crack 
propagation as a function of K t /K u determined at both 25 and 1316 °C is depicted in 
figure 8. The overall trend of the angle of crack propagation was that 9 C initially 
decreased quickly with increasing K { /K u and then decreased monotonically at K t /K u > 3. 
No significant difference in 9 C between 25 and 1316 °C was observed. It is noted that the 
scatter in 9 C , for example, 9 C = 45 to 79° at K t /K u = 9 (pure mode II), seemed more 
appreciable in the coating material than in dense advanced ceramics, 22 which exhibited a 
typical scatter of about 15°. 

Analysis of fracture surfaces revealed no explicit difference in the fracture 
morphology of the pure mode I, mixed-mode, or pure mode II specimens. However, it 
should be mentioned that unlike dense ceramics the coating material has provided a 
challenging subject on fractography because of its unique open and microcracked platelet 
microstructure. The origin or nature of strength-controlling flaws in the coating material 
has rarely been identified with any traditional means of fractography (optical and/or 
SEM) that would be otherwise sufficient for the case of dense ceramics. 12 An improved 
means of fractography is required to better distinguish fracture morphologies in response 
to different modes of fracture of the coating material. 
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3.2. Consideration of Mixed-Mode Fracture Criteria 


3.2.1. Considered Criteria 


As shown from the crack propagation angles in figure 8, mixed-mode fracture does 
not take place in the same plane as the original crack. This noncoplanar crack 
propagation requires different fracture criteria, typically a combination of the driving 
forces K\ and K n together with K\ c and Ku c . The mixed-mode fracture criteria for an 
isotropic, homogeneous material can be broadly categorized in four areas and will be 
briefly described then compared with experimental results. 

3.2.1a. Maximum Principal Stress Criterion 24 This criterion assumes that a crack at 
instability propagates in the direction normal to the maximum principal stress direction. 
The instability condition is given by the following relationship: 

Kj c - Ki cos 3 — + 3AT cos" — sin — = 0 (10) 

lc 1 2 2 2 


The angle of crack propagation 9 C is a function of K { IK n and is given by 


0 C = 2 tan 


-l 


II 

4 AT 


- + I 


kKuj 


-ll/2l 


( 11 ) 


In case of pure mode II, the following relations can be obtained from Eqs. (10) and (11): 

0 C = 70.53°, = 0.8662 (12) 

K lc 

3.2.1b. Minimum Strain Energy Density Criterion 25 This criterion assumes that crack 
propagation would take place in the direction along which the strain energy density is 
minimum. This assumption gives rise to the following mixed-mode fracture criterion 

Klc ~7T, —\- a n K i +2 a l2 K^n +a 2 2 Kf{\ = 0 (13) 

A K — i) 


where 

flj 1 = (1 + cos 0 C )(K - cos 0 C ) 

«|2 = sin 0 C (2 cos 0 C - K + 1) 

a 22 — (k + 1)(1 - cos0 c ) + (1 + cos0 c )(3cos0 c - 1) 

with 


k = 3-4v for plain strain 
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k = (3-v)/( 1 +v) for plain stress 


where V is Poisson’s ratio. The angle of crack propagation, 9 C , is given by 


(K - 1) sin 


0 C —2 tan 


-i A. 

^TT 


\ 

/ 

- 2 sin 

2 

/ 

V 


9 C - tan 


-1 K, 


K, 


it y 


- sin20 c = 0 (14) 


In case of pure mode II, the following relations can be obtained from Eqs. (13) and (14): 


0,, = cos 



Kii c _ 

2(k - 1)" 

l 6 / 

K lc 

a 22 


(15) 


3.2.1c. Maximum Energy Release Rate Criterion. This criterion assumes that a crack 
with an infinitesimally small kink at an arbitrary direction propagates in the direction 
along the maximum strain energy release rate. 26, 29 Many investigators keeps that the 
strain energy release rate, G, can be expressed in terms of K\ and K n as follows: 

G = (’ M /yf + c \2^\K\\ + c 2iKu (16) 


where cn, c\ 2 , and c 22 are coefficients. Hyashi and Nemeat-Nasser 28 determined 
numerically the coefficients and the direction of maximum G as a function of K\IKn. 
Palaniswamy and Knauss 27 proposed an approximated mixed-mode fracture criterion on 
the maximum energy release rate as follows (that was also used by Suresh et al. 15 ) 


K, 3 
- + — 


K 


Ic 


f K ^ 

A 

\ K \c) 


= 1 


(17) 


where for pure mode II the approximation leads to 


0„ = 77.4°, — 2^- = 0.817 

A 


(18) 


The fracture criterion and the crack propagation angle on the maximum G were used 
based on Eqs. (16) and (17) in this work. 

3. 2. Id. Empirical Criterion , 30 ’ 31 The empirical fracture criterion can be expressed by the 
following generalized, simplified form 


f 


V 


Ki 


K 


Ic ) 


p 


+ 


f v \ q 

A II 


Y^ic 


= 1 


(19) 


where p and q are parameters (typically p,q = 1-2) to be determined from experimental 
data through a curve fitting, and 
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( 20 ) 



with 0<y<2, an experimentally determined value. Despite the lack of its theoretical frame 
work in some cases, the empirical fracture criterion can still provide a much simplified 
representation (as database, too) of mixed-mode behavior of a material particularly when 
other (preceding) fracture criteria do not describe reasonably well the actual material’s 
behavior under mixed mode. The maximum G criterion (Eq. 15) or the coplanar crack 
propagation criterion, 32 in fact, is one particular type of the empirical fracture criterion in 
its form. 

3.2.2. Comparison with Experiment 

Figure 9 compares the afore-reviewed four different mixed-mode fracture criteria 
with the experimental data (that were presented in fig. 5), where K u normalized with 
respect to K ic , K n /Ki c , was plotted as a function of K\/K\ c for both 25 and 1316 °C data. 
The prediction for the minimum strain energy density criterion was made using a value of 
V = 0.2. Neither the maximum principal stress nor the minimum strain energy density nor 
the maximum energy release rate criterion was in reasonable agreement with the 


1.4 


1.2 


Asymmetric test 


Experiment 
A 25 °C 
O 1316 °C 



Maximum energy 
release rate 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


K\'K\c 


Figure 9. — Plots of K^/K^ as a function of K|/K| c for mixed-mode 
fracture criteria and experimental data determined for plasma- 
sprayed Zr0 2 -8 wt% Y 2 0 3 thermal barrier coatings, v is Poisson's 
ratio. 
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experimental data. Only the empirical criterion was in good agreement with both ambient 
and elevated-temperature data, resulting in the following parameters: 


p = 1; q = 2; y = 0.65 for 25 °C 
p = 1; q = 1.3; y = 0.65 for 1316 °C 


(21) 


Previous studies on dense ceramics such as silicon nitrides and alumina 22 and ceria-doped 
ziconia 23 showed that the minimum strain energy density criterion resulted in best 
agreement among the mixed-mode criteria mentioned above, as that was the case in 
another study of alumina by Suresh et al. as well. 15 Hence, there exists a remarkable 
contrast in mixed-mode behavior between the dense and the porous TBC ceramics. The 
unique nature of the coating material — possessing porosity, microcracks, and platelet 
microstructure — might have been responsible for the ability to distinguish its unique 
mixed-mode fracture pattern, as this is not the case for dense ceramics. As also seen from 
figure 9, the three mixed-mode fracture criteria — the maximum principal stress, 
minimum strain energy density, and maximum energy release rate criteria — predict the 
values ofK llc /K lc = 0.8 to 1.1. The poor agreement with these three fracture criteria was, 
therefore, due to the lower values of K Uc /K lc that were around 0.65 for the coating 
material at both temperatures. Hence, the value of K a JK lc is the parameter with the 
largest influence on the degree of agreement with the prediction of the fracture pattern. In 
case of K\\JKi c ~ 1, statistically reasonable agreement would be found in any of the major 
mixed-mode fracture criteria mentioned. The dense ceramics examined in previous 
studies 22 ' 23 all exhibited a value close to KnJKi c ~ 1. Comparison of KuJK lc among other 
brittle materials will be presented and discussed below. 

The data on crack propagation angle shown in figure 8 were compared with the 
predictions made with the three mixed-mode fracture criteria, and the results are 
presented in figure 10. The minimum strain energy density criterion overall seemed to 
yield a better prediction compared with the other two criteria. However, because of 
somewhat significant scatter in 9 C as well as little difference in prediction between the 
criteria, it is difficult to state which criterion gives the best agreement with the 
experimental data. Note that a significant discrepancy was found particularly in the 
region close to pure mode II, e.g., Ki/Ku < 1. This discrepancy predominant in the region 
of Ki/Kn < 1 was also observed in the previous studies with silicon nitrides and alumina 22 
and other grade alumina by Suresh et al. 15 By contrast, this discrepancy seemingly unique 
in ceramics was not found in amorphous glassy polymer PMMA (polymethylmetha- 
crylate) mixed-mode tested in either asymmetric four-point flexure 23 or uniaxial tension 25 
configuration. A good agreement of propagation angle with the fracture criteria (as well 
as little scatter in 9 C ) in a wide range of K t /K u close to pure mode II was also found in 
steel tested in asymmetric four-point flexure at -196 °C to induce brittle fracture. 33 

As previously mentioned, the overall locus of mixed-mode fracture (in K n versus K h 
or Kn/Ki c versus Aj/W Ic relation) is mainly controlled by the value of KuJKi c ; hence, the 
ratio of mode II fracture toughness to mode I fracture toughness is a very important 
measure of a material in terms of its response to mixed-mode fracture. With this in mind, 
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90 



Figure 1 0 . — Plots of crack propagation angle as a function of K\/K\\ 
for mixed-mode fracture criteria and experimental data determined 
for plasma-sprayed Zr02-8 wt% Y2O3 thermal barrier coatings. 


a summary of data on K a JKi c for various brittle materials compiled by Munz and Fett 31 as 
well as data on dense ceramics from previous studies 22 ' 23 was used and compared with 
the TBC data obtained in this work, as shown in figure 1 1 . The values of K U JK\ C for the 
dense ceramics two silicon nitrides (one with fine and one, elongated grains), coarse 
grained alumina, 22 and ceria-doped tetragonal zirconia polycrystal (Ce-TZP), 23 all tested 
in asymmetric four-point flexure using naturally sharp precracks, are invariably around 
KnJKi e = 1 and are much greater than those of the TBC {K n JKi c ~ 0.65). Except for this 
comparison, a large variation of K n JK\ c varying from 0.6 to 2.0 is noted depending on 
materials, test specimens, notch or initial crack preparations, and even investigators, etc. 
As a result of this significant variation, a reasonably unified trend on the value of K Uc /Ki c 
cannot be made for the materials provided. However, it should be noted that the value of 
KiiJKi e for the TBC evaluated in this work is still at the lower end of the data pool. One 
thing to note is that test technique must be technically sound, particularly in the areas of 
specimen configuration and precrack preparation. A diametral compression technique, for 
example, would be in danger to overestimate greater A IIc as a result of frictional 
constraint between the two crack faces by the existence of compressive stress therein. 22 
Use of sharp precracks using an appropriate method such as single edge precrack beam 
(SEPB) method 14 or SEVNB method as used in this work is also crucially important to 
obtain as accurate fracture toughness as possible. 


NASA/TM— 2003-2 12 1 85 


16 



1.4 

1 1 

1 

1 1 1 


1.2 

- 


- 


1.0 



^ ▲ 

o 

as 


▲ 

▲ 

1 

0.8 

“ 


Dense ceramics 


0.6 

□ D 
1 1 

Previous studies [22,23] 


0.4 

This study 




~ 


~ 


0.2 

- 


- 


0.0 

(3) I I 

i 

1 1 1 


TBC TBC 

> 

N> 

o 

GO 

Si3N4 Si3N4 Zr02 



2.2 

25 °C 1 31 6 ° 

C 

(HP) (Sintered) 






2.0 

A 


From Nunz and Fett [31] 


1.8 

- 


- 



A 


V 


1.6 

- 






¥ 


1.4 

- 


- 

o 

as 

1.2 

- 


8 s 

6 

1.0 

A 

□ 

□ 

v o * 


0.8 

A 

□ 

o 


0.6 

o & 




t 


o 


0.4 

(This study) 




0.2 

- 


- 


(b) | | 

1 

Ill 


0.0 

TBC Al 2 0 3 

Si 3 N 4 

Zr0 2 Glass Graphite 


(Zr0 2 -Y 2 0 3 ) 

Figure 1 1 . — Comparisons of K\\ c /K\ c of plasma-sprayed 
Zr0 2 -8 wt% Y 2 0 3 thermal barrier coatings with other 
advanced monolithic ceramics and brittle materials at 
ambient temperature, (a) Comparison with previous 
studies 22 ’ 23 , (b) Comparison with data compiled by 
Munz and Fett 31 . 


3.3. Effect of Directionality 

To determine a possible effect of the material’s directionality on fracture toughness, 
mode I fracture toughness testing for the coating material was performed by the DCB 
method in the direction perpendicular to the plasma-sprayed direction, as mentioned 
in the Experimental section. The value of fracture toughness was found to be K ic = 
1.04±0.05 MPa Vm . This value of AT Ic was very similar to that of K\ c = 1.15±0.07 MPa Vm 
determined in the direction parallel to plasma-spraying direction. Hence, the directionality of 
the coating material in response to mode I fracture toughness was insignificant. The previous 
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studies showed that even the tensile strength of the coatings with different vintages was 
almost independent whether the plasma-spraying direction which was either 
perpendicular or parallel. Therefore, a notion that mechanical properties of TBC would 
be strongly dependent on the direction of plasma spraying may not be supported and 
generalized, based on the fracture toughness testing in this study and the strength testing 
in previous studies. 12 

3.4. Effect of Sintering 

Figure 12 shows the results of mode I fracture toughness testing for the coating 
material annealed at 1316 °C in air with annealing times ranging from t = 0 (as-sprayed) 
to 500 h. Fracture toughness increased significantly at / = 5 h, increased monotonically to 
t = 100 h, and then reached a plateau at t = 500 h with a value of K ic = 2.6±0.2 MPa Vrn , 
about a 120 percent increase from t = 0. This increase in K\ c was attributed to a sintering 
effect. 6 Figure 13 shows typical fracture surfaces of specimens tested, subjected to 
annealing for t = 0 and 100 h. The sintering effect was manifest by the evidence of 
significant grain growth for the specimen annealed for t = 100 h. It has been shown that 
sintering gives rise to an increase not only in fracture toughness but in elastic modulus, 
strength, and thermal conductivity of the coating material. ’ 4 ' 12 ' 6 Although not performed 
in this work, it is expected that K\\ c in response to sintering would follow a trend similar 
to K lc . Therefore, evaluation of Kn c as a function of annealing time at 1316 °C is an 
immediate task for future study. Since a change in mechanical properties (and physical 
properties as well) would occur inevitably upon sintering, particularly in a short period of 
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Figure 1 2. — Result of fracture toughness as a function of anneal- 
ing time for plasma-sprayed Zr0 2 -8 wt% Y 2 0 3 thermal barrier 
coatings annealed at 1316 °C in air. 
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Figure 1 3. — Comparison of fracture surfaces of fracture toughness-tested specimens of 
plasma-sprayed ZrC> 2-8 wt% Y 2 O 3 thermal barrier coatings, (a) As-sprayed, (b) Annealed 
in air at 1 31 6 °C for 1 00 h. 
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Figure 1 4. — Typical stress-strain curve for plasma-sprayed Zr0 2 -8 wt% 
Y 2 O 3 thermal barrier coatings in compression 11 . 


time, a continuous use of as-sprayed mechanical data in elevated-temperature 
applications would give an over-conservative estimate, leading to erroneous results in 
reliability and/or life predictions of components. 

3.5. Other Considerations 

Finally, it should be mentioned that the coating material has shown nonlinearity and 
hysteresis in its stress-strain relation, irrespective of loading configuration: tension, 
uniaxial or biaxial flexure, or compression. 12 A typical result showing such nonlinearity 
and hysteresis exhibited in compression by the as-sprayed Zr 02 -Y 2 03 coatings is 
presented in figure 14. This nonlinearity in constitutive relation was due to the coatings’ 
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unique microstructure (porosity, microcracks, and platelets) and diminished upon 
annealing to a degree depending on annealing time and temperature. A question arises as 
to whether linear elastic fracture mechanics (LEFM) approach could be applicable to this 
nonlinear, elastic, as-sprayed coating material. The degree of nonlinearity was more 
significant when the level of applied stress was increased. In most cases, however, 
fracture force employed in fracture toughness testing was relatively low (with about one 
tenth of the peak force in figure 14), and in this case the nonlinearity became negligibly 
small. As a result, an LEFM approach (a continuum approach as well) would be justified 
at least in fracture toughness testing for the as-sprayed coating material. 


4. CONCLUSIONS 

A full range of mixed-mode fracture behavior covering mode I and mode II for 
plasma-sprayed Zr 02-8 wt% Y 2 O 3 thermal barrier coating was detennined in air at 25 and 
1316 °C in asymmetric four-point flexure in conjunction with the single edge v-notched 
beam (SEVNB) method. The following conclusions were made: 

1. The mode I fracture toughness was found to be A! lc = 1.15+0.07 and 

0.98±0.13 MPaVrn , respectively, at 25 and 1316 °C. The respective mode II 
fracture toughness was K llc = 0.73±0.10 and 0.65±0.04 MPaVrn . Hence, the 
coating material exhibited an insignificant difference in either A) c or K nc 
between 25 and 1316 °C, whereas it exhibited noticeable difference between K\ c 
and A) Ic , resulting in Kn c /Ki c = 0.65 at both temperatures. 

2. The empirical mixed-mode fracture criterion was in best agreement with the 
coating’s mixed-mode fracture behavior among the four mixed-mode fracture 
theories considered. The angle of crack propagation was in reasonable 
agreement with the minimum strain energy density criterion. 

3. The mode I fracture toughness of the coating material at 25 °C was almost 
irrespective of the plane either parallel (typical) or perpendicular to plasma 
spraying direction so that the coating’s directionality in response to K [c was 
insensitive. 

4. Sintering at 1316 °C in air showed a significant influence on mode I fracture 
toughness, giving rise to an increase in K ic at the plateau region by 120 percent 
when the coatings were annealed for the duration of 500 h. 
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